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Abstract 

A nonlocal method to obtain discrete classical fields is presented. This technique relies on 
well-behaved matrix representations of the derivatives constructed on a non-equispaced 
lattice. The drawbacks of lattice theory like the fermion doubling or the breaking of chiral 
symmetry for the massless case, are absent in this method. 
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1 Introduction 



The non-perturbative implementation of QCD uses the Feynman path-integral approach 
to yield a discretized version of the theory on a equispaced discrete space-time lattice 
where the main properties of the quantum fields are expected to be preserved Jl|]. How- 
ever, properties such as chirality and uniqueness of lattice fields are not preserved. In 
fact, the Nielsen-Ninomiya theorem states that for a local and translationally invari- 
ant hermitian discretization of the fields it is not possible to have simultaneously chiral 
symmetry and uniqueness (no doublers). Many attempts have been addressed to circum- 
vent the restrictions of this no-go theorem (see for example ||). Among these, one can 
find the use of higher number of dimensions [f| , the construction of Dirac operators satis- 
fying the Ginsparg- Wilson relation 0, the use of different lattice spacings for the gauge 
and the fermion fields |J, or the proposal of nonlocal actions (see for example 0] and 
references there in). The locality of the Dirac operator is in general a desired property, 
however, there are examples of operators that are topologically proper but not local ||. 
In this letter we introduce another approach to the problem of discretization of fields. 
This is done by using a non-equispaced lattice, a projection of the partial derivatives 
that have been used to discretize the classical free Klein-Gordon and Dirac fields || and 
a discrete Fourier transform [Kj. The classical discrete propagators approach to their 



continuum forms when the number of nodes tends to infinity and a relationship between 
this technique and that of path-integrals can be established. Such a technique can be 
seen as a projection of the quantum algebras on a finite linear space yielding matrix repre- 
sentations for the partial derivatives that produce hermitian actions with nonlocal kinetic 
terms ||, [H|. This scheme yields Dirac operators that are free from some drawbacks 



exhibited in the standard lattice method as the fermion doubling, the chiral symmetry 
breaking for the massless case or the yielding of different dispersion relations for bosons 
and fermions. 



2 Discrete formalism 



In order to be concise, we point out in this section the main features of the technique and 
leave the proofs of our assertions in the references |], [H], 



11, 12, 13 



We begin by choosing a non-equispaced lattice constructed with the set of (fi = 
0, 1, 2, 3) zeros Xj of the Hermite polynomial Hn^) along the direction /i. This gives a 
total of iV = N NiN 2 N 3 lattice points in a non-equispaced hypercubic structure. Now, 
we construct the four N x N matrices 
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where 1^ is the identity matrix of dimension and is the skew-symmetric matrix 

fO, i = j, 

Wjk = 1 . . 

u ui * T 1 J ' ■ 

We also need 

S = S ® S z ® S 2 ® S x 

where is a diagonal matrix whose jth nonzero element is given by the value of the 

function exp(-£ 2 /2)#^(£) at £ = a#. 

Then, 

SD M S 1 (2) 

is a projection of the partial derivative in the subspace of functions U generated by 
products of the form u n (£) = exp(—^ 2 /2)H n (^) with £ = x M , /i = 0,1,2,3 and n = 
0, 1, . . . , — 1. Note that {w n (0}o°> ^ s a orthogonal basis for square integrable functions 
in one dimension. Thus, if ip(x) G U and \I/ denotes the iV x 1 vector of components 

% = l/>(Xq) = ^{x],xl,xf,X m ), 

ordered according to 

q = j + (k- l)iVi + (I - l)NiJV 2 + (m - l)NtN 2 N 3 (3) 

where j = 1,2..., A^! runs faster than A: = 1,2..., A^2, running faster than I — 1,2..., N 3 , 
running faster than m = 1,2..., N Q , we have that the vector ^ constructed with the 
values of = d^ip & t the site x q = (x^, xj, x 2 k1 xj) and ordered as in (^), is given by 

„ = SD^S- 1 ^. (4) 

This means that (H) can be used to reproduce exactly the values of the partial derivatives 
of functions belonging to U at the lattice sites. On the other side, if ip( x ) is a square 
integrable function (considered as a function of just one variable x M ) that is not contained 
completely in U, Eq. (HI) yields approximated values for the derivatives of ij){x) at the 



lattice points. The error can be estimated in very special cases [|13[ and become, as 
expected, related to the complement of ip(x) with respect to U. This behavior can be 
shown by considering a function of the form (p(x) = a(x)ip(x) with ip e U (we will take 
just one variable in the following), and discretizing its derivative at N nodes. Note that 
<p is not necessarily contained in U and in general, the vector <J>' constructed with the 
values 

if'(xj) = aix^'ixj) + a'ixjWxj), (5) 

is not given by SDS^ 1 ^, where $ denotes the vector constructed with the values <p(xj). 
Thus, we address our attention to the residual vector R = $' — SDS^ 1 ^. Because of 
@, $' can be written as ASDS^ 1 ^ + A'^f where A is the diagonal matrix containing 
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a(xj) along the main diagonal and A' the diagonal matrix containing a'(xj). Thus, the 
remaining vector giving the error becomes 

R = {ASDS- 1 - SDS^A + A')V = SAS' 1 *, 

where A is the matrix given by A = [A, D] + A', i.e., 

' a'(xj), j = k, 

A jk = < 



a(xj) — a(xk) 
? 

Xj Xj» 



By rewriting the commutator as 

DA = AD + A' - A, (6) 

we can see the way in which the derivative of the product of functions a(x)ij)(x) is handled 
by our method: just construct the diagonal matrix with nonzero elements a(xj) and apply 
the product SDAS^ 1 to the vector The error R can be determined by considering the 
action of A on an arbitrary vector V. To this end, we take first the case where a(x) is 
taken as the power x n . It is not difficult to show that in such a case, AV is given by the 
linear combination 

A V = 6 £n-l + &l£n-2 + • • • + &„- 2 £ + 6„-l , (7) 

where £j is the vector whose jth component is x l j and the coefficients are given by h\ = 
J2jLi xl jVj- The same procedure to get (|7|) and the expansion of a(x) in powers of x 
yields again (J7|) for the partial sum a n (x) = Y<k=o a n xk i but the coefficients change to 
h = J2j=i Ai{xj)Vj, where Ai(x) is a polynomial of degree I. In any case, R is always in 
the subspace spanned by (0) and determined by the higher power of x in a n (x). Now, 
let us consider the part of R given by AV, where V = S'~ 1, I / . Since ip(x) G U, the 
components of ^ are of the form tyj = exp(—x"j/2)f(xj) (where f(x) is a polynomial of 
fixed degree), the coefficients of (0) becomes 

N f(x) 

bl^^M^Jjrr-y / = 0,l,---,n-l, (8) 

j=l n N\ x j) 

where n is the higher power of a n (x). Therefore, SAS~ X projects ^ on the subspace 
spanned by the independent vectors constructed with the values of exp(— x 2 /2)x N+l ~ 1 
at the nodes. Thus, the error in the discretized derivative of a function a{x)ij}{x) not 
belonging completely to U depends on the components Ui{x) = exp(— x 2 /2)Hi(x) defined 
by the orthogonal complement of U. Now, to show the condition under which bi = 0, let 
us take the higher power of Ai(x) in (]8|). This gives a sum of the form 

£4/M 0) 

7 -=i H' N {xj) 
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By using the matrix representation of d/dx given in |14| applied to the vector formed with 
the values of the function g(x) = (x — Xj) l+1 f(x) at the nodes, one can show that if f(x) is 
a polynomial of degree at most N — l — 2, the sum given in (|J) is zero for a power of Xj less 
or equal to /. Thus, the derivative of the product a n (x)ip(x) can be represented exactly 
by our discrete scheme whenever ip(x) is restricted to be in the subspace generated by the 
basis Ui(x), I — 0, 1, • • - , N — n — 1. In the case of four variables, a(x) can be expanded 
in terms of the basis u q (x) = u m (x )Mi(x 3 )M/ c (a; 2 )Mj(x 1 ), where = exp(— £ 2 /2)if n (£) 

and (||) becomes 

D M a = aD M + a„ - A M , (10) 

where the projection matrix A M can also be given in terms of the one-dimensional matrix 
through a similar relation to ([!]). Such a matrix makes a projection on the orthogo- 
nal subspace to the one generated by the basis u q (x). Therefore, we have a convergent 
discretization scheme if is large enough. 

Let us take a look on the similarity transformation of (|j) when — > oo. First of all, 
note that (|]) implies that D M can be considered as a representation of for tranformed 
vectors $ of the form $ = S _1 \l/. Since 



the jth nonzero element of takes the asymptotic form (-l) N ^ +j 2N^\ /T[{N^ + l)/2)], 
therefore, the jk element of S^D^S' 1 differs from the jk element of only by (— iy +k . 
This means that in the limit — > oo, the transformed vector $ = S" 1 ^ differs from \1/ 
only by an alternating change of sign along the /i-direction. In other words, if we take into 
account the fact that the similarity transformation, given by S, transforms ^ into itself, 
except for an alternating change of sign along each direction in the asymptotic limit, we 
always may use D M instead (0) as representation of <9 M . 

The commuting matrices (P can be diagonalized simultaneously by the unitary and sym- 
metric matrix 

F = F* ® F 3 ® F 2 ® F u (11) 
(* means complex conjugate) where 



N u -1 



(F,) jk = £ ( 12 ) 



1=0 

and 



(iV M -l)!2^-i-i Hi(0 



Here, pj is also a root of H N ^) and it represents an eigenvalue of the discretized linear 
momentum. It is possible, on physical grounds, to construct ( |T2"D as a convergent discrete 
Fourier transform fllQf . Denoting by [•••]& the fcth column of a matrix, we have that 

[F»] k ^ (-1) N » +1 [F% (13) 
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as — > — p^. Besides F^ = F 1 , the basic properties of the four-dimensional discrete 
Fourier transform ([0]), are 

iD F = FP°, —iDjF = FP J (14) 

and 



A' 



lim F qq , = C N e~ ip ^' (15) 



where C N is the product of constants of the formQ C Nfl = 2 N ^ 2 (T[(N I1 + l)/2]) 2 /^!, 
and p q ■ x q i = Pm x m' ~ Pf x f' ~ Pk x t' ~ p] x }'- The order between indexes is given by (|]). 
In (14|) P M is given by tensor products of the form ([!]) in terms of the diagonal matrix P M 
containing the eigenvalues p^ along the main diagonal. 



3 Discrete Lagrangians 

The path-integral quantization of scalar and spinor fields uses the classical Lagrangian 
of the corresponding field. Thus, our aim in this section is to give the rules to obtain a 
discretized version of such a Lagrangian and their properties. In order to do this, let us 
take for instance the Lagrangian density for QED: 

L = 4>{ip- ey^Afi - m)V> - \{d,A u - cU M ) 2 , (16) 

where ip = ip(x) and = A^(x) are the spinor and gauge fields respectively. 
To get a discrete form of this Lagrangian we choose the four-dimensional lattice x q of the 
previous section, i.e, zeros of -ffjv M (0 as lattice points for each direction /i = 0, 1, 2, 3 
giving thus a total of N — N0N1N2N3 nodes. To maintain the relationship between the 
continuous and discrete field variables, let us denote by \I/ the spinor vector of components 
ty aq = ty a (x q ), where a = 1, 2, 3, 4 is the spinor index and q is the lattice index ordered as 
in (^). A component of this vector can be described by only one index r according to the 
order: 

r — j + (k - l)JVi + (/ - l)iVyV 2 + (m - 1)N X N 2 N 3 + (i - l^A^A^o (17) 

where the slower index is i = 1,2,3,4. Thus, ^ is a vector of dimension 4N. The gauge 
vector, denoted by A^, has the components A m = A^(x q ) ordered as in @, yielding a 
vector of dimension N. The values of the derivative d^A v at the nodes are denoted by 
(AvJg- Therefore, a discretization of fllTf ) is 

C q = M aq jZ b (*,») bq ~ *a q (ej% b A m + m5 ab )V bq - -[{A v ^) q - (A^) q ] 2 , (18) 

1 In fact, there is a factor (— iy+ k in the asymptotic form of {F^)jk that should be taken into account 
in Cat . However, this factor can also be treated along the same lines as in the previous paragraph. 
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where ^f aq = ^Iqlba- According to our procedure, the derivatives (^^)b q and (A U}tl ) q can 
be approximated by (D u ) qr ^br and (D u ) qr A ur , respectively [D^ is used instead of @]. 
Here, a,b = 1,2,3,4 and q, r = 1,---,N. Repeated indexes are summed excepting q. 
When is written in matrix form, the reference to the index q disappears and the 
equation becomes the inner product 

S = f [iy ® (D M + ieAj - ml w ]f - ^DM" - D"^) T (D^ - D^), (19) 

where .4^ is a diagonal matrix with the values A m along the main diagonal, ^ is the 
transposed conjugate vector 

* = [(7 ® l JV )*] t = *t( 7 °® ljv), 

l n is the identity matrix of dimension n and the superscript T means transpose. Because 
(19) involves the sum over all q and this means to sum over the spatial and temporal 
indexes, S is a discretization of the classical action save for a constant which gives the 
correct dimension of the action.^ Since is real and skew-symmetric, S is real. By 
making the definitions 

D M = D M + ieA» (20) 

and 

S M = -\{^A V - D"A"f(D^A v - D^), 
S can be written in a simple form 

5=*(ip-ml4tf)* + (21) 

with the understanding that the slash notation involves tensor products in the discrete 
case. 



3.1 Chiral symmetry and absence of fermion doubling 

As we mentioned it, our discretization scheme presents some desirable features to be 
considered in studies of field theories such as QCD. First of all we want to show that the 
doublers are absent in this scheme. Let us begin with the free Dirac equation, 

(ip-m)i/) = 0. (22) 

According to our method, the discretized version of this equation is 

(17" <g> D M - mU N )^ = 0. (23) 

2 Indeed, we are using dimcnsionless quantities for both the sites at the lattice (space-time and mo- 
menta) and the discrete fields. 
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This is an eigenvalue equation for m and it can be solved by using the discrete Fourier 
transform F given by (jTH) . Denoting by M the diagonal matrix containing the eigenvalues 
m q , ( ]2"B| ) becomes 

where ^ is the matrix whose gth column corresponds to the eigenvector Multiplying 
this equation by I4 (g> F^ and using Eqs. (0) we obtain the transformed equation 

(^®F t D /1 F)$ = (7 At ®P M )$ = $M, (24) 

where 

$ = (l4®F f )* (25) 

is the matrix of transformed eigenvectors and P M = g IJtv P v . This equation, the discrete 
version of ^ip = mi/j, says that M is similar to the matrix P defined by 

P= 7 ^P M . 

By applying P again to (f2"4|), one can shows that M 2 becomes diagonal and equal to 
I4 <S> m 2 , where 

m 2 = (p0) 2 _ ( pl )2 _ (p2) 2 _ (p 3 )2> 

Thus, the square of a mass eigenvalue, say m r , is given by the discrete form of the energy- 
momentum relation 

< = (A? - (p)f - (Pl) 2 - (PIT, (26) 

where the indexes are ordered according to fli~7|) . Therefore, the doublers are absent. 
Let us note that (|26D yields a dense set of points m 2 in the set of real numbers since p£ 
are zeros of the Hermite polynomial if^^). 

In the case where there is a gauge field coupled to a spinor field, the discretized Dirac 
equation is 

= ml 4ff *. 

This equation defines the Dirac operator as iYf). Due to the properties of tensor products, 
this matrix satisfies trivially the discretized condition for chiral symmetry: 

( 7 5 ® 1 N ) (if)) + (zjZ>) ( 7 5 ® 1 N ) = 0. 



3.2 Plane— wave solutions for the free Dirac field 



Let $ 9 be 



We can obtain solutions of the free equation (|23|) by anti-transforming (|2] 
the eigenvector corresponding to m q (a vector of zeros everywhere save for the gth element 
which is one). Since m q is degenerate, i.e., there are several values of giving the same 
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eigenvalue m 2 ,, the solution of ( p3|) \I/ mr corresponding to the mass eigenvalue m r is given 
by the linear combination 

V mr = E c i®l F U (27) 

where p 2 = {p^) q {j)^)q and c 5 is 1 x 4 vector. According to the asymptotic form of the 
elements of [F] q given in ([15|) , fl27|) corresponds to a plane-wave discretized solution of 
fl2"2]). Since ( P^| ) gives the same value m 2 under the interchange p^ — > — p^ and [F] g changes 
4i 



to [F'] 9 , except for a minus sign [cf. (0)] that can be included in the linear combination, 
Eq. (|27| ) splits in two sums, one corresponding to positive values of p° m and the other to 
negative values, i.e., 

v|> mr = E + ^ ® [F], + E~ b « ® [ Ft L' ( 28 ) 

where the + (— ) superscript means to sum over points lying on the upper (lower) hyper- 
boloids defined through 

P°m = ±H + (p)) 2 + (pl) 2 + (pf) 2 ] 1/2 - (29) 

To find out the explicit forms of a q and b q , let us take just one summand of each sum of 
(|28f) and define them as 

K = a « ® [ F ^ ^ = h i® [ Ft L- 



By applying (|23| ) to ty q and taking into account that D M [F] g = — ip m \F\ q and D M [F^] 9 = 
ip M [F^] q , where p m is the gth nonzero element of the diagonal matrix P M , we obtain that 

(YPm ~ m q l 4 )a q = 0, (7%g + m q l A )h q = 0, (30) 

which are the equations for Dirac spinors evaluated at p m and m q (an expected result 
since the discretization scheme presented in this letter is not involved with the spinor 
space). Therefore, the standard procedure to find the solutions of ( |3~0|) can be followed 
step-by-step yielding the expansion of the wave function of a particle of mass m r in terms 
of plane waves 



= E + [(mr/pi) e «} ® m* + e~ ik/pS.) e «m ® t Ff 

p2=mp (=1,2 p2_ m 2 (=1,2 

where and 6^ are complex numbers and 

uj = C mr [l, 0,pf/(^ + m r ),p-f k /(p° m + m r )} T , 
u 2 = C mr [0, l,pj k /{p° m + m r ), -pf/{p° m + m r )} T , 
v] = C mr [pf/{p° m + m r ),pt/{p m + m r ),l,0} T , 



9 
V 2 



Cmrkfe/ + "V), ~pf / (p„ + m r ), 0, 1]' 



Here, C mr = [(p^ + m r )/2m r ] 1//2 , pi = pj ± ip\ and j,k,l,m are the indexes of the 
components of p q satisfying (^9|) for the given eigenvalue m r . 



3.3 Gauge invariance and the covariant derivative 

Let us consider now the way in which the discretized action (|T9| ) changes under the gauge 
transformation 

which takes the discretized form 

*^(1 4 ®E)#, \^\ + A,» (31) 

where E is the diagonal matrix whose elements are given by exp[— iea(x q )]5 qr and A jAt is 
the vector whose gth component is d^a(x g ), ordered according to @. We address first 
our attention to the term ^{jJf)— ml w )\l/ of (^) which we denote by Sdm- Under this 
transformation, Sdm takes the form 

Sdm = ® E f )[i7 M ® (D^ + z'ei^ + ieA,,) - ml m ](l 4 ® E)tf, 

where A )A , is the diagonal matrix whose principal diagonal is the vector A jM . Therefore, 

5 DA/ = #[17" <g> (E f D M E + ie^ + zeA iM ) - ml 4iV ]*. 

In this case, ( |10D takes the form 

D^E = ED M - ieEA jA1 - A„ 

and this gives 

S DM = mU N )^ - ^{i'f ® E f A M )^, 

so that, if Np — > oo, the residual vectors vanish yielding invariance of Sd a/ under fl3~T1). 
To show the invariance of Sm, the term of (|2l|) depending only on the gauge field, let us 
write the vector A jA1 of (31) as A jjU = D M A + where A is the vector of values a(x q ) and 
is the residual vector. Since a(x) can be expanded in the basis u q (x) and D M is an 
exact representation of <9 M for these functions (save for a similarity transformation), 
must vanish when N —>■ oo, so that (|3lD takes the asymptotic form 



^^(1 4 ®E)^, Afj, —> + D^A. 

Because of [D M , D u ] = 0, 5m becomes invariant in such a limit and this completes the 
proof of gauge invariance of (pi]) when — > oo. 

Now, we want to show the way in which the commutator of the discretized form of the 
gauge covariant derivative, Eq. ( p0|) , is related to the electromagnetic field tensor F^. 
By applying the commutator [D M , D„] to we get 



where we have used again (|10|) . Thus, the expression for in terms of the commutator 



of the covariant derivative can be recovered when — > oo. 
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4 Final remarks 

Two facts about the discretization procedure presented in this letter should be noticed. 

1. Since it is a nonlocal method (the matrix representations of the partial derivatives 
involve values of the field at each node) based on a non-equispaced lattice and 
the finite number of nodes breaks down the translational invariance, the Nielsen- 
Ninomiya theorem does not applies in this instance. In fact, this technique yields 
the right results: the chiral symmetry and the dispersion relation are obtained 
correctly when the number of nodes N is finite; the quantum commutators, the free 
propagators, the gauge invariance and the commutator of the covariant derivative 
are recovered when N tends to infinity. The translational invariance is also recovered 
in such a limit since the nodes becomes dense in the space-time. 

2. Such a technique involves matrices with a well-behaved structure. Thus, it could 
be a good issue to obtain analytical solutions in some cases. 
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